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SUMMARY 

An approximate method for rapid determination of the pressure 
change for subsonic flow of a compressible fluid under the simultaneous 
action of heat transfer, friction, rotation, and area change is de- 
veloped. In the development of this method, the momentum equation was 
approximated and rearranged for a convenient solution employing charts. 
This report presents both the analysis involved in simplifying the mo- 
mentum equation and the charts necessary for obtaining particular solu- 
tions. The charts provide a step-by-step solution which converges to 
an exact solution as the number of steps is increased. An illustrative 
example and comparison with more rigorous numerical solutions with con- 
ditions typical for air-cooled turbine blades are included. These com- 
parisons show that the solution converges rapidly to provide good 
accuracy. 


INTRODUCTION 

The effective design of ducts to accommodate air flow requires par- 
ticular solutions of the momentum equation for determining the pressure 
changes encountered. Simplified methods which give an accuracy suffi- 
cient for engineering purposes are in demand. An approximate method 
for determining such particular solutions has been developed for the 
one -dimensional flow of a compressible fluid under the influence of heat 
transfer, friction, rotation, and area change. 

A number of studies of one -dimensional flow of a compressible fluid 
have been published. Reference 1 presents a rather complete treatment 
including a form of the momentum equation with heat transfer, friction, 
and area change which is suitable for numerical integration. As an aid 
for expediting the numerical integration, coefficients of the differ- 
entials in the momentum equation of reference 1 have been tabulated as 
"influence” coefficients in reference 2. The analyses of references 1 
and 2 have been extended in reference 3 to include rotational forces. 
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Reference 3 also provides a form of the momentun equation suitable for 
numerical integration and tabulates "influence" coefficients for ex- 
pediting the calculations. However, in all cases , the numerical inte- 
grations as employed in references 1 to 3 are tedious and time consum- 
ing. Particular solutions suitable for rapid determination of pressure 
changes for one -dimensional flow with both arbitrary and specific heat- 
input distributions have been presented in references 4 to 7. However, 
these solutions include only heat transfer and friction. In addition, 
reference 8 developed an approximate solution for constant flow area in- 
cluding rotation as well as heat transfer and friction. The approxima- 
tions of reference 8 were based upon experimental data for cooled tur- 
bine blades tested under static conditions. Although the time required 
for a particular solution of the momentum equation and subsequent eval- 
uation of the pressure change was appreciably reduced when the approxi- 
mate solution of reference 8 was utilized, the time required for a large 
number of particular solutions, as frequently required for design pro- 
cedures, remains very lengthy. In addition, reference 8 does not account 
for flow-area changes frequently encountered in actual practice. 

Since a method for the rapid determination of the pressure change 
through ducts with heat transfer, friction, rotation, and area change 
is in demand, the momentum equation has been approximated and rearranged 
for a convenient solution employing charts. The purpose of the present 
report is to present the analysis involved in simplifying the momentum 
equation so that a graphical solution is possible and to present the 
charts necessary for obtaining particular solutions for subsonic flow. 
Heat transfer, friction, rotation, and area changes are all included in 
the solutions. The charts provide a step-by-step solution which con- 
verges to an exact solution as the number of steps is increased and is 
limited only by the graphical accuracy. Comparisons are also made with 
the solutions determined from the more exact numerical integration of 
reference 3 for typical examples of air-cooled turbine blades. 


ANALYSIS 

The general momentum equation for one -dimensional flow is approxi- 
mated and rearranged to form the basis for the construction of charts 
suitable for the determination of the Mach number change with flow under 
the simultaneous action of heat transfer, friction, rotation, and area 
change. A relation between the Mach numbers and pressures is finally 
determined such that charts may be employed to evaluate the pressure 
changes in a duct once the Mach number changes are known. The energy 
changes may be arbitrarily specified in terms of the total-temperature 
changes and may be evaluated by a procedure such as that given in ref- 
erence 3. 
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Evaluation of Mach Numbers 

The differential form of the general momentum equation for one- 
dimensional flow is (ref. 3) 


cn 

CO 

^4 

O 


dF 


dP - F^pg dx + pV dV 4- = 0 


(i) 


where x increases in the flow direction (all symbols are defined in 
appendix A) . The differential drag force can he expressed by the 
equation 


dF d = f 




- 4f 
dx = =— 


PV" 


A dx 


( 2 ) 


and the body forces per pound of weight flow, as resulting from rota- 
tion about a central axis, can be expressed by the equation 


2 


cd r 


g 


(3) 


If the differential momentum change is replaced by the equivalent ex- 
pression 


pV dV = d(pV 2 ) - V d(pv) (4) 

and equations (2) to (4) are substituted into equation (l), the momentum 
equation, as integrated between stations x and x4Ax, becomes 



'(P y )x+Ax 

(pv). 


V d(pv) - 


£Z! dx = o 


D h 2 


( 5 ) 


The first two terms of equation (5) can be integrated directly. How- 
ever, the last three integrals of equation (5) cannot be evaluated 
directly and therefore are approximated in order to obtain forms which 
can be conveniently solved. The third integral, which exists only in 
the case of area change, and the fifth integral, which expresses the 
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drag forces, may "be conveniently approximated "by using the integrand as 
the average of the values at the end points. That is, 


(pV) x +Ax 

V d(pV) 

(PV) X 


(V. 


+ w 


[(pV) 


x-tAx 



and 


Ox+Ax 






«£ 

pV 2 

dx ■> ’ 

4f x-+Ax 

fp y2 ) + 4fx 

(Pjl) 

L 

2 

2 

_Ph , x+Ax 

\ 2 Ac4Ax ®h,x 

\ 2 /xj 


( 6 ) 


( 7 ) 


In order to facilitate a convenient solution, the fourth integral, 
which expresses the rotational forces in equation (5), will he approxi- 
mated by either 


x+Ax 


co 2 rp dx * arr^^Ax 


(8a) 


or 


'x-fAx 


CD 2 rp dx - ^ 2 r x4Ax p x+Ax Ax 


(8b) 


depending upon whether the solution is initiated with the conditions at 
the station x or the station x-f-Ax. Since the integrand of either 
equation (8a) or (8b) is approximated by the. value at a single end point, 
the approximation is in general less accurate than those given by equa- 
tions (6) and (7). Equations (6) to (8) become 1 identities in the limit 
as Ax approaches zero. 


Because the density and velocity variations from stations x to 
x+Ax are approximated by equations (6), (7), and (8), the energy equa- 
tion must be employed only to determine the total energy change. 

When the approximations given by equations (6), (7), and (8a) are 
substituted into equation (5) and the first two terms of equation (5) 
are integrated, the approximate momentum equation becomes 
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F x+Ax F x + ( p ^ ^x+Ax “ ( p ^ ): 


V x + V x-fAx 


^(pV) X 4Ax " ( p ^)xJ 


2 A Ax 


4f 


x+Ax [ pV 


r 2 \ 


4f\ 


- 


®h,x+Ax \ 2 ) x4Ax ®h,x V 2 /x_ 


pV‘ 


= 0 


(9a) 


However, if equation (8b) is used instead of equation (8a), the approxi- 
mate momentum equation becomes 


Px-tAx ~ F x + (pV^)x4Ax " (pV^)x - 


V x + V X4AX 


OT 


Ax 

4f x4Ax | 


I + 4f * 

2 

_°h,x-tAx 

V 2 ) 

x+Ax ^h,x 


[(e v W - (pv>x] - 

= 0 (9b) 


pV_ 

2 


If the variables of equations (9a) and (9b) are changed by employing 
the continuity equation, the relation between the total and static tem- 
peratures, and the equation of state, the two forms of the approximate 
momentum equation can be reduced to the useful relations (appendix B) 


f ( m x+Ax> \+Ax) - 


/ ^x+Ax \ 


\ ^ J L 


F(M X ,A X ) + 


O 

orr Ax 


G(M X ) 


(10a) 


and 


F(Mx,4c) = 


M 


T x+Ax 

\ 


ySc+Ax / 


F (^x+Ax> Ac+Ax^ 


orr 


c+Ax^ 


T^f- <>(*W 

■ L x+Ax 


where 


F(M,A) = 


M 


A 7 ^ t 1 - £ H 

1 + H_±. M 2 ) 

2 J 


(10b) 


(11) 


G(M) = 


(i ♦ ifi « 2 ) l/2 

gRM 


( 12 ) 
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(13) 



(14) 


The total-temperature ratio in equations (10) may be evaluated from 
a solution of the energy equation such as presented in reference 3, 


Evaluation of Pressures 

By applying the continuity equation and the equation of state, the 
weight flow can be expressed as 


Since the relation between the total and static temperature as given 
by the energy equation is 



F MA 


(15) 


or 



1/2 


(16) 



(17) 


T 


2 


equation (16) can be rewritten in the form 



(18) 


Using P'/P = (T'/t) F ^ F to express equation (18) in terms of the 


total pressure gives 



( 19 ) 


3270 


NACA TN 3150 


7 


Equations (18) and (19) are plotted in figure 1 for r = 1.40 with 
w/\/T^/PA and w/^/P'A as the ordinates and M as the abscissa. 
From the Mach number the value of either w^/t^/P'A or w^/T^/PA can 
be determined from figure 1. Then from the given conditions and area, 
P or P 1 can be computed. 


CALCULATION PROCEDURE 

Equations (lOa) and (lOb) represent the final forms of the momentum 
equation to be used in solving for the Mach number at either station x 
or x-hAx. If the conditions are known at station x, the Mach number at 
station x+Ax is determined from equation (lOa); while if the conditions 
are known at station x+Ax, the Mach number at station x is determined 
from equation (lOb). Increasing values of x, corresponding to a posi- 
tive Ax, must be in the flow direction. A given duct can be divided 
into an arbitrary number of increments in Ax (not necessarily of the 
same length), and the Mach number change through that duct can then be 
computed by progressing stepwise from the end where all conditions are 
known. The size of each increment should depend upon the accuracy de- 
sired. Subsequent examples will indicate the convergence of the solu- 
tion as Ax is decreased. 

Both F(M,A) and G(M) have been plotted in figure 2 for air with 
X = 1.40 to expedite the solution. Since x has a secondary effect 
on the Mach number change and the variation of x encountered in air 
flow is small, the consideration of other values of x is unwarranted. 

The procedure for determining the Mach number change and thus the 
pressure change with air flowing in a duct is to use equation (lOa) or 
(lOb) and figures 1 and 2. For example, the procedure for determining 
the Mach number at station x from the conditions at station x-fAx is 
to first evaluate (eq. (14)). From the value of the Mach number 

at station x+Ax and 4 x+Axj the values of F(M X 4 Ax> 4c+Ax) and 
G(M x4 ^) are read from figure 2. Then, with the specified values of 

a) ^ r x4-Ax^/^'x-hAx anc ^ V^x+Ax/^x (^x/^x+Ax) > is calculated 

from equation (lOb). From FfM^,^) and A x (eq. (13)), the Mach num- 
ber at station x is read from figure 2. Once the Mach numbers are 
known, the pressure change can be determined from figure 1. To illus- 
trate this procedure, a typical sample calculation for an air-cooled 
turbine blade including heat transfer, friction, rotation, and area 
change has been carried out in table I with the passage divided into 
two equal parts. Double columns (column 9 and columns 20 to 27) are 
included in the sample calculation setup in table I so that the calcu- 
lations could proceed from either the x or x+Ax station. For each 
of these double columns, the left column is applicable if the solution 
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is initiated at station x-hAx, while the right column is applicable if 
the solution is initiated at station x. The first 11 columns of table 
I give the conditions which must be prescribed. Although the friction 
coefficient f is assumed to be the same at both the x and x-hAx 
stations, it may be allowed to vary if greater accuracy seems warranted. 
The remaining columns of table I give the computations required to de- 
termine the Mach number at either station x or x+Ax. The sample cal- 
culation shown is for the determination of the Mach number at station 
x from conditions at station x-fAx. 


RESULTS AM) DISCUSSION 

In order to evaluate the accuracy of the solution presented herein, 
calculations were made for conditions typical of those occurring for 
air-cooled turbine blades (l) by use of the charts and (2) by the method 
of numerical integration presented in reference 3. The results of these 
calculations employing equation (10b) with 1, 2, and 4 steps with incre- 
ments of equal length are tabulated in table II and compared with the 
results obtained from the more detailed numerical methods of reference 
3. Three cases are illustrated in table II: (l) heat transfer and 

friction; (il) heat transfer, friction, and rotation; (ill) heat transfer ) 
friction, rotation, and area change. For the first two cases, the air 
temperature was assumed to vary linearly from the inlet to the exit. 
However, for the third case, the air temperature was varied exponentially 
to correspond with flow inside a duct heated by crossflow and with con- 
stant inside and outside surface heat-transfer coefficients. The total 
pressures used for evaluating the percentage differences between those 
obtained by the two solutions were computed from figure 1 and the known 
conditions and areas. 

It is shown in table II that the solutions converge rapidly. In 
addition, the comparison between the two solutions shows that good agree- 
ment is obtained by using only two steps. Even for the third case, the 
computed total pressures agreed within 4 percent for only two steps. In 
this case the disagreement was essentially eliminated by using four steps. 
Solutions were also made with eight steps, but for the conditions used 
in the examples of table II, the improvement over the use of four steps 
was probably not within the accuracy of the charts. The approximations 
given by equations (6) to (8) are least accurate at the higher flow Mach 
numbers where slight changes in the total pressure cause large changes 
in the Mach number (fig. l) and, therefore, large changes in the veloc- 
ity and state conditions. For this reason, the Mach numbers used in the 
examples of table II are intentionally high. It should be pointed out 
that the accuracy of the approximations employed in equations (10) are 
dependent upon Mach number changes, and not the passage length, when ro- 
tational forces are neglected. Thus, the accuracies shown by the results 
of table II are indicative of the accuracies which would be obtained for 
longer passages with similar Mach number changes due to heat transfer, 
friction, and area change alone. 
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In many cases the air temper at ure varies exponentially with x as 
in case III of table II. In order to check the effect on the Mach num- 
ber change of using the more easily obtainable linear air temperature 
distribution instead of an exponential variation, the examples presented 
in table II were recomputed as follows: (l) Cases I and II, which were 

calculated for table II with a linear air temperature distribution, were 
determined for an exponential distribution, and (2) case III was com- 
puted for a linear variation instead of the exponential air temperature 
distribution used for table II. The results of these calculations showed 
that the differences between the two variations were negligible and prob- 
ably of the same order of magnitude as the accuracy of the charts. For 
this reason, it appears that although the air temperature might vary ex- 
ponentially with x, for convenience a linear variation may be assumed 
for evaluating the temperatures needed for the incremental subdivisions. 


CONCLUDING REMARKS 

A method is developed for the rapid determination of the Mach num- 
ber change and pressure change for the flow of a compressible fluid in 
a duct with heat transfer, friction, rotation, and area change. The 
method uses a step-by-step integration which, when applied to specific 
examples typical for air-cooled turbine blades, converged rapidly and 
provided good accuracy with, at most, only four steps. In addition, for 
these specific examples, a negligible difference occurred between using 
a linear or an exponential temperature variation along the duct length. 
This result indicates that the more convenient linear temperature vari- 
ation may be used for evaluating the temperatures for the incremental 
subdivisions if the actual temperature varies exponentially, as fre- 
quently occurs in practice. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 14, 1954 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report : 

A flow area, sq ft 

b duct length or span, ft 

hydraulic diameter, 4A /z, ft 

F^ generalized body force, lb/lb of coolant, ex ? r /g (eq. (3)) 

F, drag force, lb 


F(M,A) 


( 1 + I r M2 ) 


M 


^ [l - ^ (A - 4)](eq. (ll)) 


friction coefficient 


G(M) 

g 

l 

M 

P 

P' 

R 
r 
S 
T 
T ' 



r-i 

2 



gHM 


(eq. 


( 12 )) 


standard acceleration due to gravity, 32.2 ft/sec^ 

wetted perimeter, ft 

Mach number relative to passage 

static pressure, lb/sq ft abs 

total pressure with respect to passage, lb/sq ft 

gas constant, 53.3 ft-lb/(lb ) (°R) 

radius , ft 

surface area, sq ft 

static temperature, °R 

total temperature with respect to passage, °R 
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w 

X 


X 


X-hAx 


p 

2 


velocity relative to passage, ft/sec 
weight flow, lb/sec 

distance measured from passage inlet in flow direction, ft 
ratio of specific heats, 1.40 


4f x Ax 


D 


- 2 


h,x 


x 


Ax + Ax 


1 (eq. (13)) 


4f x+Ax Ax 


D 


- 2 


h,x+Ax 


( Ax+Ax \ 

\Ax " / 


(eq. (14)) 


mass density, slugs/cu ft 


2 A 
0 ) rAx 


0 ) angular velocity, radians/sec 

Subscripts : 

e exit of passage 

i inlet of passage 

x distance measured from passage inlet in flow direction, ft 
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APPENDIX B 


DERIVATION OF APPROXIMATE MOMENTUM EQUATION 
Upon collecting terms and employing the relations 

, (PV) X (Vy.) (pv 2 ) x (pVW 

( ° V, x4Ax VT/ * ( pV >x<Ax \T/ 2 ( D v) x 


and 


(pv). 


^X+Ax^ / ( p ^)x+Ax /^x+Ax\ ( p ^ ^x+Ax ( p ^x 


x \ 2 / 


= (pv) 


(pv), 




2 TpVJ 


X+Ax 


equation (9a) becomes 


^x+Ax ^x 


(pV 2 ) 


X4Ax 


- 2-2 


(pV) x ^x+Ax^* 


IpvT 


x+Ax ^h , x+Ax 


(pvO 


2\ r 


X 


„ „ < pV W , "x^ 

- 2 - 2 ”(pv) v ■ + S 


h,x _ 


- 0 ) r p Ax = 0 
xx 


and equation (9b) becomes 


P x+Ax p x 


(pv 2 ) 


x+Ax 


- 2-2 


( pV ) x 4f x+Ax^ 


IpvT 


x+Ax D h, x-tAx 


(PV 2 ), 


„ „ ( pV )x4Ax , 4f x Ax 

■ “7^7 + 


- “ r x-tAx p x4Ax Ax = 0 


If the continuity equation is employed to obtain the relation 

w = (pVA) x = (pVA), 


'x+Ax 


or 


) x ^x+Ax 

TpvT 


'x+Ax 


(20a) 


(20b) 


(21) 
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and. if the simplifying notations 




^,3 


Ax+Ax 


A = 

X 


x+Ax °h,x+Ax \ ^x 

2 _ <x> 2 rAx 


4 ^x 4 Ax^ x g (VfAx ^ 


(13) 

(14) 


are used, equation (20a) becomes 


P x4Ax " P x " 


(p v2 ) X 4Ax , . ,, (f> v2) x ,, ,, 

4 ^ ^x+Ax - + 4 (A x “ 4) - 2 >c T x p x 


= 0 


(22a) 


and equation (20b) becomes 


P - P 
x+Ax x 




, T 1 n =0 
x+Ax x+Axx+Ax 


(22b) 


A more convenient form of equations (22a) and (22b) can be obtained if 
the variables are changed by employing the equation of state, the energy 
equation, and the continuity equation. From the equation of state, 
pv 2 can be expressed as 


pv 2 = JL v 2 = rPM 2 (23) 

gRT 

In addition, the static pressure can be eliminated by using the con- 
tinuity equation and the equation of state as follows: 


or 


w - gpAV 



w RT _ w I T 'R [JT_ 1 

a v _ a /\|rg iJt- M 


(24) 


( 25 ) 
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where the total- to static -temperature ratio can be eliminated by using 
the energy equation in the form 


11 

T 


1 + 


y-i 

2 




(17) 


to give the expression 


w |Tji 
^ M (l 


Finally^ application of the continuity equation gives 


(26) 


T ,p _ vT'p _ 1 

gpAV gA V 


(27) 


T 'p 


_w_ I t 1 WL I 

Ag 'JrgR N T M 


or, by incorporating equation (17) in equation (28), 


(28) 


T’p 


_w_ [ T^ 
Ag l\Jr gR 




M 


(29) 


If equations (23) , (26 ), and (29) are substituted into equations (22a) 
and (22b) and w /y/R/yg is factored out, the two forms of the approxi- 
mate momentum equation become 


4 


T 


x+Ax 


a x+Ax 


M (l + nil pj2 ") 
x+Ax. y 2 x-fAx ' 


V 2 


rM x+Ax 

i - — - 4 ) 


Ax 


«x ( 4 + ^ 


rM; 


i --r(Ac - 4 ) 


Vt 7 (l + ^ M 2 ) 1 ^ 2 

_ 2,. J_4 A 2 = o 


A x gR 


M„ 


(30a) 


and 


/ 
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X4AX 


■^x-tAx 


*x 


“xl 1 + ^ 4) 


«x^(l + *§i 4*ta 


rM x 

1 - (Ax -4) 


172 


rM 

1 - (A x ^ - 4) 


72 


‘x+Ax 


V T x+Ax ( 1 + T 2 ^^Ax> 

■^X+Ax®^ ^X4 Ax 


Upon simplifying, rearranging, and employing the definition of 2, 
tion (30a) becomes 


F ( M x4Ax> At+Ax) “ 


1 


T x 

[ ^x+Ax ) 

^x+Ax 

\ *x >1 


o^r^jAx 


and equation (30b) becomes 


F ^ M X-Wix' A X4iX^ 


T X4AX 


8(*W> 


where 


f(m,a) = 


M 


.( 




172 


['-4 


(A-4) 


G(M) = 


r-i 2V/ 2 

gRM 


I l/2 

= 0 

(30b) 

equa- 

(10a) 

(10b) 

(ID 


( 12 ) 
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TABLE I. - SAMPLE CALCULATION AND SETUP 


Step 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


(10) 

(n) 

(12) 

(13) 

(14) 

(15) 


b, 

ft 

Ax, 

ft 

Ax> 
sq ft 

^x-hAx > 
sq ft 

D h,x> 

ft 

^hjX-fAx' 

ft 

°R 

T x+Ax> 

°R 

r x4Ax> 

ft 

r x> 

ft 

rad/ sec 

f 

A*/ A x+Ax> 
(3)/(4) 

^X-hAx 
A * 

1/(12) 

rp i ) 

x+Ax 

(7)/(8) 

4fAx 

D h,x’ 

4(H)(2) 

(5) 

1 

0.3 

0.15 

0.000409 

0.000194 

0.00668 

0.0044 

970 

1160 

1.417 


796 

0.0065 

2.108 

0.4744 

0.8362 

0.5838 

2 

0.3 

0.15 

0.000590 

0.000409 

0.00749 

0.00668 

763 

970 

1.267 


796 

0.0065 

1.443 

0.6930 

0.7866 

0.5207 


Step 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

4fAx 

^X-tAx’ 

4(H)(2) 

(6) 

(15)-2[(12)-1] 

^x-tAx' 

-(l6)-2[(l3)-ll 

^Sc+Ax 1 T x 

Ax TiW 

(13) (14) 1 / 2 

c ° 2r x4Ax Ax 

m i > 

1 x+Ax 

(10) 2 (9)(2) 

(8) 

o^r^Ax 

m t > 

X x 

(10) Z (9)(2) 

(7) 

Mx+Ax 

«x 

G (*w) 

xlO 3 

G(Mx) 

xlO 3 

1 

0.8864 

-1.632 

0.1648 

0.4337 

116.1 


0.629 


0.963 


2 

0.5838 

-0.3653 

0.0302 

0.6146 

124.2 


0.2105 


2.80 



Step 

(23) 

(24) 

(25) 

(26) 

(27) 

^^x+Ax' Ac+Ax^ 



o^r^Ax 

£ F (*Ac4Ax> Ac+Ax ^ “ 
(23)-(24) 

[ F (Mx>Ac) + 

co2 r .Ax H 

T r ° (Mx) > 

(23)+(24) 

r(25)i 

LTWj x 

F ^ M x4Ax^ A x+Ax^ > 
[(25)d9H x+Ax 

Mx 

^x+Ax 

T« 

x+Ax 

(GtMx^) 
(20) (22) 

T x * 
(G(^)), 

(20) (22) 

1 

2.345 


0.1118 


2.233 


5.148 


0.2105 


2 

5.027 


0.348 


4.679 


7.613 


0.1350 
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TABLE II. - COMPARISON OF RESULTS FROM NUMERICAL AND 
GRAPHICAL SOLUTIONS 


Case 

T! 

X 

A. 

X 

Exit 
Mach 
number , 

Me 

Inlet 

Mach 

number 

Mi 

from 
ref . 3 

Inlet Mach 
number M^ 
from charts 

(P 1 ) - ( P * ) 

' i'ref ^ i' 

chart 

T 7 

e 

A 

e 

(P 1 ) 
''i'ref 

Trial 

Trial 

1 

2 

4 

1 

2 

4 

a I 

0.7595 

1.0 

0.833 

0.419 

0.406 

0.416 

0.418 

0.0310 

0.007 

0.003 

b II 

.6864 

1.0 

.800 

.505 

.472 

.498 

.501 

.0510 

.009 

.006 

C III 

.658 

3.04 

.629 

.141 

.1232 

.135 

.140 

.134 

.042 

.009 


^eat transfer and friction. 

b Heat transfer, friction, and rotation. 

q 

Heat transfer, friction, rotation, and area change. 


3270 


Weight -flow parameter , — pt^” ° r 
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Mach number. 


M 



Figure 1. - Evaluation of weight -flow parameter from Mach number 
(ratio of specific heats, 1.40). (A 1577 - ”by 17-in. print of 
this fig. is attached.) 
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Mach number 

(a) Mach number range, 0.045 to 0 # 080. 

Figure 2. - Evaluation of F(M,A) and G(M) (ratio of specific 
heats, 1.40). (A 12- by 20- in. print of this fig. is 
attached . ) 
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o 

E 

CO 

to 



(Id) Mach number range, 0.075 to 0.25. 


Figure 2. - Continued. Evaluation of F(M,A) and G(M) (ratio of 
specific heats, 1.40). (A 142- “by 20-in. print of this fig. is 
attached. ) 





P(M,A) 


22 


NACA TN 3150 



.35 .40 .45 .50 

Mach number 

(c) Mach number range, 0.20 to 1.0. 

Figure 2. - Concluded. Evaluation of F(M,A) and G(M) (ratio of 
specific heats, 1.40). (A 14^- by 20- in. print of this fig. is 
attached . ) 
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Mach number, M 


Figure 1. 


Evaluation of weight-flow parameter from Mach number (ratio of specific heats. 


1 . 40 ). 
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17.0 


20.8 


20.6 


20.4 


20.2 


20.0 


19.8 


19.6 


19.4 


19.2 


19.0 


18.8 


18.6 


18.4 


18.2 


18.0 


17.8 


17.6 


17.4 


17.2 


17.0 



.05 

Mach number 


(a) Mach number range, 0.045 to 0.080. 


Figure 2. - Evaluation of F(M,A) and G(M) (ratio of specific heats, 1.40) 
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12.8 


12.6 


12.4 


12.2 
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11.8 


11.6 


11.4 
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10.8 


10.6 


10.4 


10.2 
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(b) Mach number range, 0.075 to 0.25. 


Figure 2. - Continued. Evaluation of F(M,A) and G(M) (ratio of specific heats, 1.40). 
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(c) Mach number range, 0.20 to 1.0. 


Figure 2 


Concluded. Evaluation of F(M, A) and G(M) (ratio of specific heats, 1.40) 
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